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iiBBTRAGie 

A forward scatter technique can he efiectively 

used to probe troposphere for various pararaeters such 

as cross path -^jind velocity, isotropy/anisotropy of 

t'urbulence, layer heights, wave number dependence>-- 

exponent 'p* of turbulence spectruin and structure 
2 

function with moderate variations. This can be 
achieved over extended comon volume in space* 

This thesis is the theoretical study of power 
received in troposcatter communication links based on 
layered-atmosphere model. The common voluEie is 
systematically moved through space by swinging receiver 
and transmitter antenna beams along as well as across 
the great circle path, 

Por this study, the atmosphere is assumed to 
possess layered-structure. The tropospheric conditions 
are simulated by considering various models of single 
t-orbulent layers of different vertical widths at 
different heights, multiple layers, with different 
suggested values of exponent ’p' and system response 
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obtained,. Acciirate scatter signals are evaluated by 
detailed estimation of common volume duly considering 
complex forward scatter geometry along tlie great 
circle path and at a horizontal distance from it, 
.Although, specular reflection does exist and is super- 
imposed on the scatter signals, it has not been 
included in the simulated model. 
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CHAPTER I 


IHTROPUCTIOH 

1 Historical BFickgroimd ; 

irom tlie first ever recorded troposcatter commu- 
nication during 1932-33 by Marconi over a distance of 
168 miles using 550 iHz, to tiie numerous present day 
models y attempts liave ‘been afoot to explain precisely all 
the observed facts in a troposcatter system. 

The earlier models put forward to explain the 
propagation of short radio waves beyond ralUlio horizon 
were based on 

(a) Diffraction arround smooth earth 

(b) .Reflection from elevated layers caused 
bjr abrupt changes in the temperature and 
hujiiidity of the earth's atmosphere and 

(c) Duct propagation ■ 

All these models were not sufficient to satisfac- 
torily explain propagation of radio waves between 50 IHIz 
to 5 Oliz; over distances upto 1000 Kms. This led to the 
belief of some sort of scatter mechanism. This was 
further confirmed by the icresence of fields at all times, 
though subject to rapid and severe fading. 

The first theory of scatter mechanism was ’oased on txM 
theory of scattering from irregularities of dielectric 
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constant of tropospiiere caused by turbulence (Booker and 
G'ordon, 1949) • Later, Yillars and V/eiskoff (1955) 
put forward different physical causes for refractive index 
fluctuations. Obiikhoff’s statistical theory of tm-bulence 
was applied to explain the problems of radio wave 
scattering (bilverman, 1956). 

Although no single theory exists which explains 
all the observed facts of a troposcatter system, Tatarski's 
theory of wave propag'ation through turbulent medium 
(Tatar ski, 1961) provides by far the most plausible 
explanation. 

1.2 4tmo snher i c Hi cx o s true tur e ; ' , 

Ihe propagation of very high frequencies through 
is 

lower atmosphere^contr oiled by the vai^iations of the 
refractive index 'n'. i’he refractive index is in turn a 
function of temperature, pressure and humidity, 

Por frequencies less than 1 00 G-Hz^ the value of 
refractivity If ( = (n-1 )x1 0^) is given by 

N = 77.6 I + 375000 -V (1.1) 

correct upto 0 . 5^0 for 

Pressinre P between 200 and 1100 mb 
I'emperature f be tween 240^ and 31 0°K 
Partial pressure of water vapour less 

^ p ^ 30 nb. 
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It has recently been shomi by the use of high 
resolution microwa've and accoustic radars^ that there exists 
a wide variety of turbulent eddies having thicknesses 
ranging from 1 m to -a few tens of meters. Large amplitude 
.Kelvin-Helmholtz (K-H) waves have also been reported in 
the middle and upper troposphere. These are attributed 
to the dynamic instabilitjr of sheared stratified flow 
when itichardson number (Hi) ^1/4 (B irkemeiyer , 197^)* 

The forward sca'iter reflectivity is dependent on the scale 

of eddies given by 

X 

L^ = sin ©/2 (1.2) 

where X, is system wave length and © is the scattering angle. 

•1^' nay be cal.led as "System-filtering scale". 

At normal frequencies of operation (tipto a few GHz) the 

system is dependent on large scale turbulence structures. 

The micro scale turbulence which is also present in the 
doe s 

layersj/ not contribute significantly to the received signal. 

But for shorter wave length, the scatter would be due 
to the smaller scale eddies. These occur more at higher 
altitudes. The fading, in such a case xfould be faster 
since the scatterers at higher altitudes are sub jec ted to higher 
wird velocities. This also produces greater Doppler shifts. 
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Exact knowledge of re fr activity obtained from 
in-situ measiirements is essential for obtaining tropo- 
spberic microstruct-ore for correct prediction. Unfor- 
tunately the microwave refractrometers, for this purpose, 
are not available in our coimtry. Iloxfever for 2,1 GHz 
EMPUit-lTAlHIlAli link, set up by I.I.T, KAI^PUE., some in-si tu 
measurements have been carried out (Hathur, Hehra and 
iiaman, I974) by flying the refractrome ter over the 

proi^agation path. I'liese measurements have revealed that 

is 

the average value of scale of turbulence/of tens of meters 

— ^ 

and mean squared refractivit;?- fluctuations (/\n ) about 

-11 

10 for heights upto about 2300 meters. 

^ *5 Problem Defining ; 

A troposcatter link at 2,1 GHz has been esta- 
blished between KaHPIIR. and hAlhHAl ( a distance of 
330 Kms) and a lot of data has been logged over this link. 
Initial simulation for this link was developed by using 
the concept of atmospheric structures consisting of la^^ered 
scatterers (Gupta, 1 S 73 ) • Later a study of troposcatter 
signal characteristics, based on the random layered scatter 
model, assviming that the power received is iDroportional to 
(-1 1/3) th power of the scattering angle was carried 
out,(Bassi, 1974). 
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Several troposcatter pa^a^ieters such, as cross 
path wind velocity, scattering layer height, isotropy/ 
anisotropy of turbulence, wave numbei- dependence, exponent 
(p) of the spectrum of turbulence, including specular reflection^ 
contribution to the signal level can be determined by 
means of scatter experiments involving measurement ' of 
spectrum of received signal, average doppler shift and 
total power received as a function of the antenna beam off-set . 

i’or .vanpur-l^ainital link, for various assumed values 
of the above referred parameters, doppler cross wind 
relations have already been calculated and plotted (mathur 
and Mehra ,1976). 

Avery large number of troposcatter experiments 
can be used as remote probing tools for tropospheric 
parameters (Gj.essing, 1969). These fall within 7 main 
Categories. 

1) Experiments in beam swinging' in horizontal 
and vertical planes. 

2) Comparison of time-average power on several 

radio f re quencieo. 

3) Measurements of vertical correlation distance 
of a field strength. 

4) Measurement of horizontal correlation distance 
of the field strength. 
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5) Ileasurement of antenna to medinjji coupling 
loss. 

6) lime-delay experiments (Rake-techniques). 

7) Band width, measurements (frequency sweep 
experiments). 

Beam swinging experiments can also he, effectively 
used to determine the moderate variations of structure 
function (Guessing, 1962). 

In the Case of various studies conducted for 
iCanpur-hainital link, the received power (f^^) has been 
estimated based on RBb 101 report. According to this 
method, the value of the path loss is given as 

jj = 30 log f, ... -20 lo.p d +R(©d)+Aa ('l«3) 

liip rirlZ 

where Aa = absorption loss and 

I'(Od) = Scatter loss, '©' being in radians and 'd' 
in Ems . 

The scatter loss is given by the relation 

F(©d)= K(R ) + 0.34 ©d + 3© log (©d) 
s 

where IT is the surface refractivity and 

; ' . S, ■ . ■ 

K(1T ) = 139.5 + 0.06 h - 2.4 x 10*"^ ijf (1.4) 

s s s 

Thus the scatter loss as per RBS method is ' 

basically a function of scatter angle '©' and separation 
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■between the ti-ansmitter and receiver. It is thus based 
on rough approximations. Ihrther approximations used for 
calculation of the common volume in order to determine 
the received power leave a lot of room for improvement. 

Attempts have thus been made, in this thesis 
to evaluate the total power received on Kanpur -hai nit al 
link calculated on the basis of the layered atmosphere 
model with detailed analysis for more exact calculation 
of common volume. I'he factors which influence greatly 
the power received, the influence of which have been 
considered for received power calculations are 

1 ) Height and thickness of the sca/L'tering layer. 

2) Presence of more than one scattering layer 
within the common volume. 

3) Beam swinging in different directions (influence 
of antenna pointing) . 

4) Sensitivity of the received powei’ to the power 
index 'P' of refractivity turbulence spectrum. 
Phis is achieved by considering the average 
value as p = 1 1 /3 with extreme values as p=2 
and 5. 
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"* Break up of the work ; 

This work lias been divided into variOLis cha^pters. 

In cha,pter II, forward scatter geometry and various 
coordinate systems suitable for various t;?-pes of scanning 
techniques are given in brief. Common volume environments with 
ray path considerations, geometry of the coimaon volume 
and its evaluation based on Lammer's method are discussed 
in chapter III. In chapter IV a brief description of scatter 
•theory along with scatter-power evaluation and evaluation 
of scatter- integral are discussed. I'he efiects of 
specular-reflection on the received signal level are included 
in chapter V . In the last chapter, discussions, conclusions 
and recoiffinendations on future development if any are 


included 



CHAPTER II 

P0R¥ARD scatter GEOHSTRY AR'D coordinate SYSTEtlS 

2 'I Porwaxd scatter geometry ; 

Tlie basic forward-scatter geometry indicating 
great circle plane geometry of the troposcattei* path is 
given in Pig. 2.1, 


In Pig. 2.1, R,,. and T-- are 'the receiver and transmitter 

iL 

locations respectively. P is the scatter point which is 
assumed to be the center of the common volume. The centre 
is defined as the inter-section point of the antenna 
bor esights. Elevation angle are and C5{ with 

andQ(g.gi indicating the azimuthal angles. © is the 
scattering angle which is defined as the angle between 
the transmitted signal ray and scattered signal ray 
towards the receiver at point P. 


The extent of scattering is determined by the 

magnitude of flu.ctuationsof tropospheric refr activity at 

~ 2' 

the scattering point ( /a R ) . The forward scattered 
power Pg^ is given as (A-t las etal, 1968). 


2 [ . ,, 

■^T X 1 (0) 

TT^.3 J r2 ..,-2 

UTT) yi^-^T 



( 2.1 ) 




1 1 


It is thus seen that the rocoived pox^er depends on 

Transmitted pox-jer P ^, More transmitted 
po-werj more xfould be the received power. 

But there are physical limitations deciding 
optimum value of P,,,. 


( b ) T he transmitter and receiver antenna fiain 
functions G'r.* Ih-e gain functions are 

: 1 K 

dependent on the si-ze and shape of the 
antennas and also on the system wavelength. 
However antenna size cannot be increased 
too much since it results in aperture-to- 
me dixim coupling loss. 


( c ) The distance of scatter point from th e 
transmitter (Lrji) and receivex' (ij^) • At 


larger values of L,, 


the scattered 


power received will be appreciably reduced. 


('^ ) Scattering angle 'Q' and pox'^rer inde sc ' p ' . 

nvcn small increase in the value of © 
reduces Pt^ considerable the values of the 
powei- index *p' is generally accepted to 
be (11/3) (Tatarski, 1361 ; bjessing, 1962), 
Hoxfever, the extreme values of 2 and 5 are 
also considered. 


(e) bcatterin^ cross section CX/yQ) . It is 

defined as the power density of the scattered 
wave in the direction of receiver, per unit 
solid angle, per unit scatter volume per unit 
power density of the incident wave on the 
scatter volume. 
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QoEi-'^on VO lime . It is in turn dependent on 0 and 
antenna bores iglit inclination with the hori- 
zontal and the antenna beam width. 


Power- received is related throug,h (Tv'C©) to the 
three dimensional power spectrum of the refractive index 
field 0^ (E ) by the relation 


TTir^ 

(P-iQ) = 0^ (K ) 


( 2 . 2 ) 


Considering Kolmogorov’s theory of locally homogeneous 
turbulence and. ^Obulnhov ' s theory of passive admixture in 
the turbulent flow, the valne of 0^ is given as . 


0^(S)=O.O33 = .033 (2h sin0/2)"‘'^/^ 

(2.3) 

K is called the scattering vector. It is defined as 

the 

vector bisecting the angle made by/transmitter 1'-,^, receiver 

^ .A. 

R and the scatter point P. Geometry showing scattering 
JL 


vector is given in Pig. 2.2. 

Ihe magnitude of the scattering vector is ■ 

given as 

_ _JlIL sin (0/2) (2.4) 

' ' >- 
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2 . 2 n. . yo lume s c a.nnin & : 

i'iiere are five- basic scan types available for 
effective mapping of the common volime through space 
(Lammers, 1970). These, are. 

(V) Vertical Scan . The common volurae is scanned 
in hei^t above a fired position on the surface of the 
earth. If P (center of the common volume) is defined by 
the co-ordinates Tp and height li above the surface 
of the earth, then this scanning technique corresponds 
to keeping Zp, Yp constant. This does not ch.ange the 
direction of K , though- its magnitude varies. 

(2) Horizontal layer scan along path . The 

coijiiuon volume is scanned in a plane pai'allel to the great 

circle plane at a lired height. It coi responds to constant 

Y._^ and H, but changing Z . 
r ■ ■ P . 

(3) iiorizontal scan trans ver-se to the •path . This 
is essentially the same as (2) above except 'iiiat it is 

in a i3lane perpendicular’ to great circle plane. It 
corresponds to keeping }Ip and H const, ant and varying X,. 

(4) Constant scatter angle scan . In a fixed 
plane containing path terminals, by varying the height above 
the earth's surface and its position on the surface, the 
com on volume is swung in the space in such a way that the 
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scattering angle is kept constant, ihis corresponds to 
keeping and constant, I'he sc,atter position 

movements in the space do not change K, It can he used 
to determine si^atial inhomogeneities. 

(5) Circular scan tran s verse to the oath . It is 
same as (4) hut is in a plane perpendicular to the great 
circle plane. In this scan the scattering vector K 
changes direction with no change in the magnitude. 

Out of these scan types, those at serial ITo. 

(1 )> (2) and (5) are used in order to completely map the 
Gonmion volume and for determination of the tropospheric 
structures. 

2,3 Coordinate System s - Reference and Local s 

Consider Pi^. 2,3* I'he great circle plane is 
defined in the usual manner hy the transmitter receiver 
it™ and C the center of the earth, fhe following assumptions 

jL' 

are made for the purposes of analysis. 

(a) A corrected radius of the earth ( =4/3 Actual 

radius R) along the great circle path is used in 
order to compensate for the average decrease 
in the refractive indez value with increase in 
height, Ihus earth is assumed to be a sphere 
of 4/3 R radius. The values of correction factor 



16 


k considering the refractive gradient over 


Delhi region averaged over a period of 12 
months comes to However, for compu- 

5*20 4 

tation purpose k = recommended by 

CGIR for standard atmosphere has been used< 


(b) The transmitter and receiver antennas are 
assumed to be located on the surface of the 
earth. 

Let T^ (X^, Z^) and 

local co-ordinate systems at transmitter and receiver 


respectively. Then from Fig. 2.3 
2D = path length 

Tj = (-Dj 0;0 ) with X,Y,Z as the reference 

■o _ / ^ co-ordinates system 

R-- - V Dj 0; 0 ; 


The mid path height of earth's surface above 


jL-J jlane would be 




2 2 

pi 1 


1 /2 


(-1-2) B 


Assuming D S-» if can bo simplified as 

.2 


H. 


D‘ 


2(-|-R) 


(2.5) 


Let a point P (assumed center of common volume) be defined 
with coordinates (Xp, Ip, Zp). The height of earth * s 
surface above 2-1 plane with reference to point P is 
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ti 


. 4 . 1 4 . n2 . 2 2 

= XL - |(— iL-xi^) + + Yp 


with. 4.^ H- y,? XL, it reduces to 

r r :> 

H' = f ii - (|ft-H ) I 1 + -i-t-iL 

3 3 o .) 4 2 

2(^ R-H^) 


1/2 


or 


2 2 
H 

° 2(j R-H^) 


A 

But Hq ^ Pl, which caXL be neglected in the denominator 


2." + y 


H ' h:: K, 


2(f R ) 


( 2 . 6 ) 


H’ represents the radial distance of earth’s 
surface from Z-Y iJlane . But because Z^, ^ . 

it Can be assumed to be the surface height. 

Any arbitrary scatter position in an actual linh 
can be referenced by its surface dist.ance from the mid 
path pointy along and across the great circle plane and 
its height above xhe surface of che eartn,xi. this corres- 
ponds to xL., Z_, (= xi' + li) in (xl,Y,Z) co-ordinate system. 

© 


and of (+ 


x^^angular rotation of (- ) for y^, Z,^) 

^ ) for Y , is required for reference 


2 ' *■'* 'Tl* R“ 

of central coordinate system (1,Y,Z) to the local coordinate 


system (Big. 2.3) • This rotation yields 
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0 © 
cos(-“ ) + i Bin ) 


T = Y 


© 

2r|i = 2 cos (-^ ) - X sin ) 


© 

{ 

2 


and 


© 


= X cos (“^ ) Z sin (“^ ) 
Y = Y 

© © 

= Z cos (-^^) + X sin <-^2. ) 


(2.7) 


Substitution of Z in terms of II, and R in (2.7) gives. 

the local coordinates in terms of reference co-ordinates as 

© © 

Xrj, = X cos (“~ ) + (| R+ Z~I-I^) sin ) 


I,j = I 

0 

Z,j, = (4 R Z"H ) cosC—p^ )_x sin (• 


© 


'i ^ 'P 

— J- tx. 


and 


© 0 
Xj = X oos a. + Z - H^) sin (-j^) 

= Y 

A ® ^ 4 

Z^'= xt +Z-xl^) cos (— )+X sin(-^2 ) _ | 


R 


2 • 4 Elevation and Azimuthal angles ; 


( 2 . 8 ) 


Erom the local coordinates of the transmitter and 
receiver locations (Eig. 2.1 ) for the given scatter geometry, 
the elevation and azimuthal angles are 





HIT 


t-^C 


7^ 


'iransLiitter azimuthal angle = iirctan (” 

u 

'iransinit ter elevation angle = Arc tan 

z... 


. “”5™— T72 


o( 


HR 


ck. 


T£L 


Pueceiver azimuthal angle = Ajrctan 


Receiver elevation angle 


Y.. 

il 


u. 


^R 


Arc tan 

■■ Z-, 
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Oii^TEii III 

COi#iOK YOWm DlTIjiONiiSMI'S 

5 • 1 - S-aY Path considerations with Diffraction effects ; 

Por evaluating any troposcatter link design, the 
detailed knowledge of the atmospheric stmicturcs and its 
effects on the signal passing through it is very essential* 
The signal ray which passes through the atmosphere suffers 
any or all of the following 

a) Scattering 

h) Reflection 
c) xiefraction 

Besides these, the phenomenon of diffraction affects the 

signal level considera'bl 3 ?' under certa.in conditions, rhe 

path 

diffraction in a tr anshorizon2 such as troposcatter link 
occurs in two ways. 

i) Jjif fraction around a smooth spherical 
earth's surface. 

ii) Ifnifo edge diffraction” 

The diffraction over a smooth homogeneous earth 
may he determined hy the first term of the Vander Pol- 
Bremmer series (Bremmer, 1949). Certain graphical methods 
are also known for the determination of signal attenuation 
in the diffracting regions for hoth horizontal and vertical 
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polarization an! ior any combination of surface conduct- 
ivity^ permittivity^ , frequency o± operation '±*, 
distance between terminal points 2D and antenna iieigiits. 
However the diffraction effects are pronounced only 
within the distances upto which the diffracted signal is 
comparable with the scattered signal, fhe distance 
beyond the radio horizon at which the scattered and 


diffracted signals are equal is given by 

r~ — r1 /3 

4 = 65 r-%^ n Km 

L J 

For Kanpur -r'ainital link operating at f 


(Ohisholm, 1 956) • 
(3.1 ) 

2.1 GHz. 


d = 23.5 Xm • 

iiius we c.an say that for our link, for distances greater 

' ' ,t ' 

than 23.5 Kms beyond the radio horizon, the scatterea field 
will doraina,te and hence diffraction effects would be 
neglected. 


In case or knife ed^e (defined as sharply 
defined obstacles opaque to and placed in the way of 
radio waves) the diffra-ction effects are significant only 
if the coaraon volume is in very close vicinity of such an 
obstacles. In such a case Fresnel-Zone clearance from, 
the obstacle top is to be considered (Dalukhanov, 1971). 
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i'or ii^axipur-lainital link, there is no such 
obstacles on the way and hence diiir action is not 
possible. 

^ ^ kav lath considerations and Goriuaon Yolu iD. e beometry ; 

ihe propagation of radio waves throug'h the 
atrao sphere may encoimter monotonicall;v decreasing refra- 
ctivity and/or" some abrupt changes in the refractive index 
with the increa.sing height. In the first case normal 
bending of rays due to changes in refractivity can be 
appro3i;imat ed by the use of modified effective radius of 
earth. The bending oi rays can also be calculated by 
simple application of bnel3. *s law. But in the second case, 
if there is sufficient abrupt change in the variation of 
the refractivG index with height (and some times with 
distance as well) stratification of the atmosphere is said 
to take place at that height. The stratifications or 
layers reflect a substantial amount of radio energy at 
the elevations and frequencies of interest. The reflection 
co-efficient of the layers decreases with increase in the 
grazing angle or frequencies. 

In the presence of layer structures, the signal 
strength is expected to be larger than for a well mixed 
troposphere. The results expected are as shovm in Big. 5»1 
(Benglor, 1971). 




FI6 3 1 RELAT/v£ 3ICA/AL SyREf^t^THS 
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The signal strengths have also been correlated 
■with the occuTahce and nat'ure of layers throiagh 'k' 
diagrat^Sj •!£;' being the effective radixjs of the earth 
factor (jPig. 3.2). The k'-diagrans are the k‘- pro file 
with height. It is soon that the value of k. increases 
to about 7.3 fron its normal values of 4/3 at a height 
of about 1 Em above sea level. It remains unaltered 
upto about 1 .5 Em height when it again comes to its 
norraal value. This sudden increase in k value provides 
necessary condition for super -refractive layer (when 
dN/dh lies between -100 N to -157N, N being r efractivity) 
(Mitra, 1975). 

The comcaon volui?.o of the Kanpur-hainital link 
is located in a hei^t range of 1.8 Km to 4.6 Ki'a. There- 
fore the recommended value of k (= ^) is used. 

iilthoiagh there is not much of a correlation 
between the field strength and k alone^f or values of 
k^. 1 .6, improvement is obtained by correlating field 
strength with a factor k'*^. which includes the height and 
thickness of the layers (Marquart, 1970). It is defined 
as 

n d .b 

Tooo ' “ 1 00 

^ m m 


k* = k - a 


(3.2) 
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•where 

*a* = a parameter representing height influence 

h = a Parameter representing the layer thiclmess 
influence 

h = actual height (of layer) 
d = actual thickness (of layer) 

hj^ = the height in which the layers have naximun 
efficiency. 

Further, scattered power is a function of 0, 

0 *“P 

in that it depends on (sin ^ ) , 'p' being power 

index. The value of +p' may lie between 2 and 5. Ihe 
value generally accepted and used is — ^ (latarski, 1961). 
I'hus any increase in scattering angle Q will increase 
the height and decrease 'bhe magnitude of comon volume, 
fhe increase in the height of common volume tends to 
decrease the intensity of turbulence causing weak 
scattering, since ZXn , mean square fluctuations in the 
refractive index, decrease with height. Hence one has to 
compromise in keeping the beams as low as possible and 
gaurd against the obstacle blockade. 

For Kanpur-Nainital link, the axes of aatennas 
tangents 

are parallel to the£ ^ to -the earth's surface at 
the point of antenna mounting instead of pointing towards 
horizon. The geometry of the' link is given in Fig. 3.3 . 
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The heights of tronsnitter and receiver antennas 

above the mean sea levels are 1892 neters and 145«8 

meters respectively, which are very small a,s compared 

to effective radius of the earth (= ^ B. = 8500 Exa) • 

Hence these are neglected and the antennas arc assumed 
located 

to he^at the surface of the earth. The modified link 
geometry therefore would he as in Big. 3.4. 

All the points located above tangent planes 
I 2 and (Big. 3.4) have the property of being simu- 
ltaneously seen from the terminals of the link. The 
collection of these points within the troposphere 
contained in the area bounded by the t.angent planes 
1 ^ , L 2 » ^3 called scatter volume or oomiion 

volume. This is true, however for narrow -beam antennas. 
In the case of antennas having larger antenna beam 
widths, the common volime is not dependent on the 
tangent planes but on the (©)~^ power dependence. 

3.3 Link parameters for Kannur-Nainit al Mnk i 

Height of Transraitting antenna above grotind level 

= 1 4 ft or 4.3 m. 

Hei^t of Transmitting antenna above mean sea level 

= 1892 m 

Height of Receiver antenna above ground level 

= 60 ft. or 18 . 4 m 
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Heigiit of Heceiver ajatenna above neab- sea level 

= 145.8 n 

Distance between transmitter and receiver 

= 330 Kn. 

3 dB bean widths of receiver and transmitter 
antennas = =0.95 degrees = 16.6 MR 

Maximum antenna gain = 42 dB. 

3.4 Estimation of common vol\xme ; 

In the study done by Gupta(1973)> the common 
volume for Eanpur-Uainital link has been calculated 
approximately. The antenna cone tangent planes anp 
assumed to constitute rhombic common volume geometry. 

The 'parameters such as 0 , mean height of common volume 
calculated as above have further been used for study of 
troposcatter signal characteristics based on random 
layered scatter model (Bassi, 1974) and for determining 
the ’Rake* characteristics of the link (Rama Rao41975).- 

But because of the assumptions of narrow 
horizontal layers that partially fill the common volume, 
backed by the experimental observations (Lammers, 1970) 
the detailed size of the common volume is required to be 
calculated. Bor this purpose consider Figs. 3.5 and 3.6, 
These give the general orientation of the scatter systems 
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with scatter plane and forward scatter geonetry which 
are used for determination of the scatter volumo 
limits to be used in the integration. 

I 

Amid path cartesian system (Z,Y,Z) is chosen 
for basic reference such that = (D;Oj 0) and = 

(•!>; 0; 0) . A scatter plane is assumed to be hinged to 
the X-axis (Fig, 3» 5)* This plane is defined as a 
plane containing R^, and scatter point P, the center of 
the common volume , 

For the sake of convenience, a subsidiary set 
of coordinates such as (X,”^! , ■^) are chosen for the 
estimation of the common volume limits. The subsidiary 
coordinate system chosen is such that its origin 
coincides with the central reference coordinate system 
( X, I, Z) , with perpendicular to and within the 
scatter plane. The transmitter and receiver antennas 
beams are replaced by cones of equal width. The 
scatter plane determines the inclination ^ of the sub- 
sidiary coordinate systems used in the integration process, 
and ©j^ are the antenna elevation angles in the 
scatter plane and © is the scattering angle. In order 
to provide complete (X, Y,Z) coverage of the scatter 
volume center point, the"^ variations are within the 
limits as 7T 
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Consider the intersection of the cones on the 
scatter plane, as through 1^. The intersection of 
(1^ > L^) and ( 1 / 2 > determines the X-limits and the 

intersection of (L 2 > 1 ^) and ( 1 ^ , L^) determines the 
^-limits. Then with the scatter volume center position 
^p) scatter plane 


hip 

9^ = arctan ) 

Sv 

®a = arotan ( 5 ;^ ) 

The tangant planes from through would he 


( 3 . 3 ) 


TR1 = tan (Oj^ + ^/2) 

TR2 « tan (9^^ - p/2) 

TT1 = tan (9,^ + p/2) 

TT2 = tan (9,^. -/^/2) (3.4) 


The linear equations for the, planes, thus are 


Por 

S = TR 1 

(D-2) 


Por 

II 

(D-X) 


For 1- 
3 

= TT1 

(D+X) 


For 

^ = TT2 

with 

(D+X) 

(3.5) 

intersection 

of these/one 

another give 

int egration 


limit s on I.-3 plane as 

Y - - TT 1 

>^min ~ TR2 + TTT * 
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■^ax 

IR1 
" TR1 

- Ti2 
+ TT2 

*D 


-Smin 

^2 

TR2 

. Tf2 
+ TI2 

. 2.D 


^■max 

II 

1 

. TT1 
+ TT1 

. 2.D 

(3.6) 

For determining"*^ 

i 

) -limits, consider Pigs. 

(3.7) 


and (3.8). Consider an elementary volume S = 

Projecting S on X -^plane is obtained. Uow any 
plane throxogh S and perpendicular to the transmitter 
beam direction would intersect PT at Por "^^>0 the 

intersection of the plane through S and of antenna 
cone of apex angle indicated hy the semi circle 

of the radius SSq determines the^J -limits. 

If is the distance of the scattering point 
from the transmitter, then considering relative distances 

^ L,j, 


and 



(3.7) 

Also 


(S„sp2 ^ (33^)2 

(3.8) 

Further 






1 s 



tan 

L^o^^l ~ D + X 



♦ ' 

Z-^0^^1 ~ arctan ( jj^^) 



1 
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or S, = arotan (^) - «, 


T 


S S, 


Thus ^ 0^1 ~ ^T 


(D+X) " ®T^ 


Similarly 

^ ^T /^T 

Substituting (3.9) and (3*10) in (3.8) 


(3.9) 


( 3 . 10 ) 


2 2 | ^2 
9^ ~ ^T I j^T ^ 


^ p 

- 


D+Z 


(3.11) 


Replacing by half antenna beamwidth gi^ves the maximum 
distance of the transmitter antenna cone above the 
scatter plane. Thus replacing by y^/2 and substituting 
(5,7) for Lrjj in expression (3.11) gives one of the 
-limits of the integration as 




1 /2 


( 3 . 12 ) 


In a similar way for the receiver antenna, the ^-limits 
are obtained by replacing O^i with ®R- /Sj and 


substituting expression 

4 = (D - +S^ 

instead of in (3.1 1 ) as 
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V) 


IVz 


max 




\ , ^ 


( "'2" ) “■ ( 9^ 9 *^ ( ]q 3 ^^ 


- V 


1 /2 


( 5 . 13 ) 


The limits for negative ^ values are 

. = ->! 

i min ■ 


max 


and 


° 'Is, 


( 3 . 14 ) 


’max 


Thus knowing all the limits the integral can 
he evaluated using any of the standard integration 
techniques. 

The angular distance of elemental volume from the 
antenna beam axes is given as 

>)2 , ^ .71 1/2 


h = 


and 


(D+X) 


+(arctan ^ - 0^)' 


h 


r-7,2 




1 /2 


(D-X)^+S 

The antenna gain functions are, therefore given as 


(3.15) 


( 3 . 16 ) 


= G^.Exp[_-f ./S^J 

and = G^.Exp £-f .y6|] (3.17) 

where G^ and f are respectively the peak gain and antenna 

shaping factor. The incremental scatter volume is weighted by 

E n ' 

(3.18) 
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CHiiPTEB IV 

DBTBRI4IlilTI0N OB SCATTER POlf{SR 

Brief description of Scatter Theories and Scatter 
Mechanism i 

Three .distinct groups of theoretical work have 
been used to explain the scatter mechanism. These are 

(a) Theories based on tiirbulence (Booker -Gordon) 

(b) Mode Theories (Eullinsten etc) 

(c) Reflection Theories (Briis, Hog, Crawford) 

Bach theory has the backing of some reasonable experi- 
mental work. 

The elements of Booker-Gordon, Villars - 

Weiskopf's and Tatar ski's theories are explained in 

fair details in the ea^rlier works (Gupta, 19735 Bassi, 

1974 ). The theory of reflection from . the infinite 

flat layers (Briis theory) is discussed in chapter 7. 
of 

Some details/scatter theory based on the work of Guessing 
are 

(G jessing, 1 9 ©) £ discussed here. This scattering theory 
determines the fine scale properties of the refractive 
index structure. The refractive index spectrum is 
defined as 

<^ (K ) = ) K j 


( 4 . 1 ) 
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where K = K - K | K , K being wave . nmabers of 
O S 0 s 

incident ahd scattered fields and 0 is the angle 
between K and K (Pig. 4.1 )» and 

OS 

I = -~^sin (| ) (4.2) 

Assuming incident wave to be plane and that 1he 

-1/3 

linear extent of common volume ( bt \7’ '■) very small 
compared to the distances between common volume and 
the terminal stations, the scattering cross-section 
is given as (G-jessing, 1 969) 

O- (Q) = H-f- <:^(K) ' (4.3) 

where C^( K ) is the spatial power spectrum of refractive 
index, fhus in general from (4.3) the three dimension 
spectrum can be expressed as in (4.1). With isotropic 
refractive index field, the power spectrum of scattered 
wave is given by 

P(©)o< (4.4) 

There can be three different cases for 
sca'tl^©^!^ mechanisx as given below. 

Case I 

Scattering by homogeneo us isotropic turbulence : 

In the wave number (K) range corresponding to 
the inertial sub range , the one directional refractive 
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refractive index spectrum E(K) can be escpressed in 
terms of turbulence theory as 


-2/3 g-5/3 


where 


<£'e 


B(K) 1 

2 

^ = standard deviation of permittivity 


(4.5) 


/ av 


fluctuations 


1 = input scale 

27T 

( = ^ , Em being wave number for 


maximum 


E(K)^ 


E(K) ,in effect, is one dimensional power spectrum 
obtained by integration of (K) over all directions 
of scattering vector. In case of isotropic medium, 
it takes the form 

E(E) = (p (K) (4.6) 


ihis gives the refractive -index spectrum in the 
inertial sub range as (G jessing; 1964a) 

<piK) = 0,03 0^ (4.7) 

where = 5. 3<;S"et^ (4.8) 

= Tatarski's structure function coefficient. 
In (4.8), p has been replaced with (-1 1 /3) in accordance 
with Tatar ski's theory of turbulence. Substituting these 
values in (4.3) gives the scattering cross-section oK©) as 


cr-ce) = 0.03 c 


n 




(sin I ) 


>-11/3 


(4.9) 
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Case II 

Scattering from stable regioa ; 


Assuming single horizontal layer through -which 

refractive index varies systematically based on the mean 

temperature mean humidity and mean pressure, the scattered 

field 1 at a distance R from scattering element, under 
s 

the influence of incident field is given by 



E, 


s 


®s “ ^ 


A (K) 


(4.10) 


where A(E) = /fS) exp ^''"^y y ■*" ^z ’2)1 *^2 

(4.11) 

and f(r) is the space f-unction (G 3 essing,i 964 a» ^ 969) . 

But it is assumed that pe,r-mittivity is constant in X 


and Y direction and changes only in the vertical 
direction. Then f(r) can be replaced with f(z) and 
therefore 




. z 


r 


d. 


-dKy.y 



-jK .x 

e .dx 


( 4 . 12 ) 


In (4.12) the integrals with x and y arguments are the 
Fourier representation of the delta function. Thus 


_ f 

A(K)= |f(z) e ^ 
J 




(4.15) 
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Tile function /^(K) and hence E » therefore 

s 

vani^ unless "both K ahd K are zero. This means 

X y 

that for 1 to exist, E must be normal to the layer 

o 

boundary. In this case scattering cross section is 
given by the product of A(K) with its complex 
conjugate as 


cr ( 6 ) = 


2 



e 


z 


#z 


dz 


2 


Cage III 


( 4 . 14 ) 


Scattering from strong turbulent layers ; 

The intense turbulence affecting propagation, 
in this Case, is confined to the horizontal layers of 
very limited vertical extent. The layer boundaries 
will not be smooth in terms of ^ , But these layers 
are associated with large degree of local stability. 
There are therefore two factors contributing to the 
scattered power; 


(a) Scattering from mean profile as discussed 
in case II above, 

(b) Super imposed on this, the scattering 
due to random refractivity fluctuations 
caused by intense turbulence. 

The resultant scatter field is sun of the two and is 
given by 
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<r'(©)= .03 (sin + 

jV 

(4.15) 

The relative importance and hence contribution of the 


- 1 1/3 ^ tr k" 


^ -jK .2 1 
f(z) e 2 ^ * 
1 .dz 


.2 


two terms depends on 

(a) relative magnitude of r.m.s. random refra- 
ctive fluctuations and degrees to which 
refractive index varies through the layer, 

(b) thickness of the layer in relation to the 
quantity X/sin 0 and 


(c) radio wave length 4^ relation to input 
scale *1’ of turbulent field. 


If the layer thickness is comparable with scale 


arr 

K 


which contributes to scattering at a particular X 


a particular ©, then the layers have pronomiced effect. 


The power received at the receiver antenna may 
be expressed 

= P +p 
R s d 

where P is power due to scattering by turbulence and 
s 

p^ is the power due to glints or facets (the glints/facets 
are the sharp boundary variations within the common 
volxme). The refractive index modulus M(= R+1 57) gradi^t 
completely describes the value of P_. In case of narrow 
beam antenna it is given by (Mitra, 1975). 
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P = P 
s 


0 7fi n2 .dM x-U/3. ■ -11/3 

fs . (2D) 


p? 

X -p 


(4.16) 


where P^^ is free space propagation power, 2D is the 
distance "between and R^, and ^ is ahtenna "beam width 


in radians. 


In case of P^, the 


troposphere is assumed 


to "be completely mixed excepts for glints. This can 
"be characterised "by an equivalent or effective refle- 
ction area A. . The geometrical area of glints is 
assumed to "be very small as compared to first Fresnel 
zone clearance "but is large as compared to the wave 
length. In such a case P^ is given by (Mitra, 1975) 


'f s 

l.S. 


H . . dv 

e 

y 


(4.17) 


In cases where random fluctuations in the 

-5/5 

refractivity follow the inertial sub range K' law 

(meaning there-by that <P (K) is proportional to the 
-11 /3 

value of K ), the dependence of received power on 
the scattering angle is 


P(G) oA (K) . V 
But for narrow bean antennas 


( 4 . 1 8 | 
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.-1 




Vc< © 

(E)o< ©“^ 

P(6)o^ 0^5 


( 4 . 19 ) 


This is true for scattered power. In case of specular 
reflection, assuming linear decrease ot J^n over the 
height interval /Vh, the power reflection coefficient 
is given by 


|p2j _ j ^n 


. 2 / Q \ 
2 sin 


/ sin X N 

. ( — - ) 


(4.20) 


where X is defined as 


.. 27rAh . © 

X =. — sini 

h 


(i) 


The two limitting cases depending on the size of the 
layer are 

(i) Infinitely thin layer (Ab-^ro). In this case 






n 


^ . 2 / Q V 

2 sin ( '^ ) 




s-4 


which gives 

P(©)r-' ©"‘^ (4.21) 

' .f , ' ... 

' ( ii) Infinitely thick layer ( A h->co) . The 

angular spectrum in this case is given by 
2i {/ dn A, 


iri 


dh 


sin (^ ) 
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which, gives 

P(0)r-/©”^ (4.22) 

4.2 Scatter Integral and its Evaluation : 

The forward scattered power arriving at the 
receiving end is obtained by integrating over the 
common volume 7 and is given by (2.1 ), where c7~^(©) is 
the scattering cross section (also called reflectivity 
of the turbulent medium per unit volume). In one form 
of the many representations of Tatarski's theory, it is 
defined by (4.9) considering polarisation factor also, 
cr(0) takes the form of . 

<r (e) = 0.58 (sin f)'"'^^^. sin^ 

(4.25) 

where J:^is the angle between the electric field vector 

of transmitted ray at the scatter point and scattered 

2 

ray towards the receiver. Sin^ accounts for the 
polarisation effects. Expression (4.25) considers 
isotropic spectrum given by eqn (2.5). 

While considering narrow beam forward scatter 
2 2 

system, C and sin v, are taken as constant for 
integrating over the common volume (Lammer, 1972). 
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^R ® considered variable for 

tJois purpose. Substituting (4.23) in (2,1 ) ^ives the 
scatter integral » 


R^=1.91 . 10”^ 


P ^ 5/5 
^T* A • n 


. 2 

sin^ 


(j 


R 


.2.2. 


jLjlJCsin |) 


£n 1173 


dv 


V 


(4.24) 


This gives theoretical value of P,, • When it is 
compared with the expreimental observations there exists 
some descrepency. This is due to the fact that 


i) The troposphere cannot be considered turbu- 
lent throughout its extent and 

ii) significant scattering taRes place over a 
limited region only. 


Based on the study of data from a large number 
of troposcatter links, the region of significant scatter 
in the troposphere can be assumed to be a thin layer 
(Hafdy et.al, 1966 ; Crane, 1968), However there may be 

more than one layer affecting the scattering from the common 
volxime. High resolution radars used to probe the 
troposphere have suggested that a horizontal layer of 
vertical extent of 100 meters is representative of 
the actual conditions, I. i. < . i / i ^'UR 

CENtKAL UaRARY 

Aik mm ' 
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In order to solve scatter integral (4.24) the 
integration is carried out over the auxiliary coord- 
inate system (xjlj, ^), The boundary limits of bounding 
surfaces are determined as per (3.6), (3.1 2) and (3«13) 

■‘^ven after determination of boundary limits, because of 
extremely complex nature of the common volume no 
analytical solution of scatter integral (4.24) is 
available. It is therefore solved numerically. For 

this purpose a limitting surface of (X- X . ) and 

■ ° max min 

height ( S) 


max 


^in) is chosen in z - ^ plane for any 


arbitrary scatter plane inclination^. If a small 
elemental volume (dx d”^ d3) within this surface is 
taken, it will contribute to the scattering only if its 
height above the scatter plane lies within the transmitter 
and receiver antenna cones* Whenever this condition is 
met by h particular volume element, it indicates the 
extent of the scatter integral involved. Por obtaining 
received power, the above is to be weighted, with reference 

to the scatter plane , by . appropriate gain functions 

2 2 ' 0 —11 
G-jp and 1^, with © dependency as (sin ^ ) '"^ . 

Depending on the position of scatter volume and that 

of the scatter plane inclination '2^ with relation to 

great circle plane, a height above the surface of the 
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earth is attached to above. A n incremental volnmes 
within the common volume are, in this way, assembled 
together into sets of height intervals for any arbitrary 
scatter plane inclination. 

There is a change of coordinate system from 
reference to auxiliary , ^ ) in order to achieve 

integration. Por this purpose consider Pig. 4.2 which 
shows arbitrarily chosen general location of the 
scatter point P. 

It is seen that 
2p = PA + AB 

and PA = -^|cosY , AB = DO = S sin "2^ 

. Zp =Ssin‘^ -yjcosT^ (4.25) 

Similarly 

= DC - BC ( = AD) or 

X 

Yp =^cos ^ + Tj sin (4.26) 

Solving (4.25) and (4.26) for and Sgires the values 
of anxillary coordinate systems in terms of central 
reference coordinate system as 

= Yp sin'^- Zp cos^f 

5 = Yp 000^^ + Zp sin'g' 


(4.27) 
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For ascertaining the height of any scatter element 
aToove the surface of the earth consider Fig. 4.3 where 
's* shows the thickness of the scattering layer. Then 

z 

tan = — Y ~**~ 

P 

and Z=:^sin'^ - "Vicos'^ 

Substituting for sin and cos in terms of Tp and 


Zp gives 


Z = 


; /p _ Vi ^P 

+ z2)172 i (i|Tz|P^ 


or 


Z = 


s.z^ ~ ]].ip 


p 


r „2 „2 'j Tj2 

jjp ^pj 


(4.28) 


In the actual computer processing all elementary scatter 
volumes are referenced by their height Y5ih with reference 
to the center of the common volimae as 


A = Z - Z. 


(4.29) 


The approximate height of the scatter element above the 
earth's surface is 
H = Z - H, 


o 


or 


H = 


SZp -^Yp 


rv2 

) ? 


+ Zp 


2 jT 72 


(D )^-Xp - 

2 ( I R) 


(4.30) 
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Thus knowing all the parameters, the value of the 
scatter integral is computed in order to determine the 
system response. 

4 *3 Polarization effects : 

Consider Pig. 4.3;^, The angle x is between 
Lg, and B and E is normal to P^,. In case of horizontal 
polarization E has components in X and Y dire ctions only. 
In case of vertical polarization E has components in 


all directions. 

Por horizontal polarization : 

Only direction of E and is of any significance. 
It is assumed that (lammers, 1970) 


_ _zl__ E = o 
E^ “ D+X 7 Z 

then 

1 = (-1; D + Xj o ) 

= (D - X; -Y; -Z) 


• ' cos}<il = 


- Y .XD 


||(D +X)^+Y^Jj[(D - X)^+Y^+?jJ^^ 


(4.31 ) 


The factor for accounting for Ihe polarization is determined 
from (4.31 ) as 




FIG. 43A.P0LARIZAT10N EFFECT ON SCATTER SIGNAL 





F or vertical polarizatio n ; 


The direction of E is such, in this case, that it 


is resolved into and Erom treangle TP'P 


E 


Z 


I' 


(D + X) + Y' 


2 , ,. 2-1 1/2 


■] 


b^y = z 


Bjj can be further separated into its const itutent 


components E^, Ey as 


and 




Ey = 


.(D + Z) E 


XY 


i^(D + X) 

- Y . E 




XI 

2 ^ ^2l^72 


Rd +x)^ + y^ T 

^ f” 2 2~ 1 /2 

Multiplying each component with |(D + X ) J ' 

yields E as 

1=1^- (D +X) . Z; -Y.Z; (D +X)^ + Y^ 

1;^= ( - (D + X), Y,Z ) 


COS: 


- z (D 4-X) ,2 -D 


I |^Z^(E+X)^+Y^j +1^ (D+X)^+?J|P^'(D+X)^+ 
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The polarization factor is therefore 

2 2 
sin-^ = 1 - cos -yu 

The angular range of polarization is 

(f- - + Q) 

For Kanpiir-Uainit al link, the value of 0 for grazing 

incidence and fixed antenna position is 38,8 MR. 

However the beam swinging in the plane perpendicular 

to the great . circle plane results in appreciable change 

2 

in © Vg_lue. Therefore the term sin is not neglected 
while computing the received power. Farther the feed 
to the antenna dish is vertically polarized for this 
link. Therefore expression (4*32) is made use of for 
calculation of polarization factor.'- 


4,4 Results of Simulation ;, 

A system response is simulated taking into 

consideration the parameters of Kanpur-Hainital link. 

The simulation program is given at appendix 'A', The 
single-layered model consists of the horizontal layer 
of 100 M vertical extent , This layer is assumed to be 
located at the center of the common volume. The beam 
swinging in the vertical plane is achieved by considering 
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the heights of layers at 2,0 Kras, 5.0 Km ahd 8.0 Kras, 
above the surface of the earth, for each layer height, 
the variations in the power received as‘ a function of 
cross Path variations are obtained by varying the value 
of y. This provides system response for cross-path 
b e am- swinging , 

The simulation also accounts for system response 
for different layer dimensions. For this purpose a layer 
of 400 meters vertical extent is chosen at 5.0 Km 
height. 

The condition of encountering more than one 
layer within the common volume is taken care of by chosing 
two layers of 1 00 meters width each separated by a distance 
of 500 meters. The layers are assumed to lie symmetrically 
on either side of the center of the comraon vol\mie. 

The power received for various layer heights 
considered for different Y-positions is plotted in 
Fig. 4.4. Fig. 4.5 depicts variations of the length of 
common VO Ixmie ^ Z as a function of scatter position. The 
values of scattering angle © in mid path and great 
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circle plane are shown in Pig. 4.6. Por a represen- 

-tative layer of 400 meters at a height of 5.0 Kms , 

the system response -for different values of exponent 
11 

P = 2, ~ aud 5 is shown in Pig. 4.7. The system 

response for two layoisof 100 meters thickness, separated 
hy 500 meters ahd situated symmetrically on either 
side of the center of the common volume, is plotted in 
fig. 4.8. The height of the common volume center point 
is assumed to ho 5 Kin above the surface of the earth. 













Layer thickness = 100 M each 



Y (km) 


FlG.4-8 SYSTEM RESPONSE FOR MULTi LAYERS 
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CHAPTER V 

SPECULAR-RE ELECTION CONTRIBUTION TO RECEIVED POVfflR 

5 • 1 R sglew of Eriia Theory of reflection from in f inite 
flat layer ; - . _ 

% 

In terms of this theory (Eriis et.al; 1957 ) 
xmcorrelated reflections from the assumed layers in 
atmosphere caused hy relatively sharp gradients of ref- 
ractive index are assumed to be responsible for propagation 
of signal beyond radio horizon. This may not be true 

that 

entirely, in that, it has been established^the scattering 
of radio waves from turbulent changes in the refractive- 
gradients takes place accompanied Dy the specular ref- 
lection from heights vrhere mean vertical refractive 
gradients show a sudden change (Atlas et.al; 1969). But 
the specular reflection, which by its very nature takes 
place only near the great circle, is some tines found 
out to be several times the magnitude of the power received 
via scatt.erors. Hence the necessity of determining the 
specular contribution. 

The reflections from any layer depends on the 
amplitude reflection coefficient a function of 

incidence angle and layer dimensions with respect 
to Eresnel-zone., dimensions. There would be three 
different cases based, on iiie size of layer considered 
(Eig. 5.1 ). 



L arge layer s; 


In this Case the layer is assumed to he a plane, 
perfectly reflecting surface of unlimited extent, which 
is' also the basic assumptions made for scatters 
(infinitely long layer with very narrow vertical extent ), 
The power received in this case is given as 



"^T * '%?. 

_ — 

4X ) 


2 


q. 


(5.1) 


where ’q[' is called the co-efficient of reflections and 
Aqi » a,TO effective areas of transmitting and receiving 
antennas. The expression (5.1) is applicable only if 


(a) Horizontal extent of layer is much greater 

than , ^/^eing the incident /reflection 

angle for the signal, and 

(b) width of the layer is greater than » 

being separation of transiuitter and 
receiver. 

Case II 

Small layers ; 

This is the case when the layer dinxinsions are 
very small when compared to the Hrosnal zone dimensions, 
but are fairly largs as compared to the wave length. In 
that Case (jFriis, 1957) 
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Am • An n n O 

where G is width of layer h is length of layer. 


(5.2) 


Case III 


of Inter mediate size: 


In this case C is laige hut (h,ip) is assumed 
ho fairly small as compared to the Fresnel-zone dimensions; 



Now many of the layers contributing to the power 
received are not necessarily in the horizontal plahe. 

They may he oriented in any direction. The reflection, 
therefore would take place from layers located on "off" 
as well aS 'ON' side of great circle path. If 7 is the 
common volume and there are N contrihuting layers per 
unit common volume, then for intermediate layer size 


Pg, ^ T * 




2^(^ 7 


N. • N7 


(5.4) 


The only assumption made in this case is that the number 


of layers per unit volume does not alter throughout the 
common volume.’ If If is^taiken as change of refractivity 
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gradient (or the change in the gradient of dielectric 
constant at the boundary lawyers), then 



The specular ref lectivity has been used 

in formal analog with (gr‘(©),the turbulent reflectivity 
for Calculation of specular reflection power contri- 
bution. By its very definition, will have finite 
value only at the position X = 0, Y = 0, Z = some finite 
Values, lying within the common volume. Off-setting 
the antenna booms off the midpath position, the 
power actually received through reflection is reduced 
by the antenna gain product along the great circle path. 
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due to specular reflection is due to reflection from the 
great circle path (laniniersj 1970). If height of the 
reflection point from the surface of the Garth is H*, 
then the reflection point will he described by the 
coordinates 



The distance of reflection point from transmitter 
and receiver locations is 

L = + (Hq + ^ (5.8) 

Considering Fig. 4.1 and substituting the values of 
expression (5.7) gives 



Substituting expression (5.9) in the expressions (3*4) and 
(3.7) gives the distance of center of common volume from 
transmitter and receiver as 
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= ( D + Xp) + Yp +(H^+H - 


2 - 11/2 

2(|r) 

J 


(5.10) 


. n. . ^,".,=.1...,- S;!. 

^ “ 2(iR) 

Then fron Pig. 5.2 


and so 


PP* = ] X? + +(H* -H+ 

I -t r o 4 -D \ 


2(f R) 


3/Sp 


2 2 

+ Lp - PP* 

2 L T L„, 


(5.11) 


( 5 . 12 ) 


(5.13) 


- PP* 


and cos^^ = — 2 l7lT( (5.14) 

Expressions (5.13) and (5.14) give tho values of^,^ ahd 
which, suhstituting in expression (3.18 ) yiold the 
ant onna gain functions product G-^ G^. This antennas - 
gain product when .suhctitutod in basic received po-roiv 
scattGr'intCi_,ral (t-.24) giver? out tho power received 


due to the specular reflection. 
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chapter VI 

DISCUSSIOH and COHGIUSIOH 

Discussion on results pHtained ; 

The model simulated to obtain the system 
response consists of assumed layer widths typically 
less than the vertical extent of the common volume. 

The received power is calculated with respect to three 
different layer heights 2,5 and 8 Kms above the 
surface of the earth. The layer height conicides with 
the centre of the common volume. Besides these, the 
received signal for the grazing incidence, which is 

the condition of the actual Kanpur -Uainital link, and 
height 1 . 6 Km above the surface of the earth (common 
volume at 3.2 Km above the x-y plane) is also computed. 

The discussion on the results obtained can be 
broadly devided into three categories, namely, (a) 
changes in scai'bering angle for different common volume 
positions in space (b) changes in linear extent of common 
volume vs the position in space and (c) effect of beam 
swinging on the received power. These are discussed 
separately. 
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Change in Scattering angle Q ; 

The beam swinging along the great circle path, 

keeping cross path position and height of common volnme 

above the surface of the earth unaltered, has very 

little effect on the scattering angle. Tor beam 

swinging upto 20 Ems the change in the scattering angle, 

is of the order of 0.1 degree only. This is shown in 

Tig 4 • 6. However with variations in cross path position 

0 changes considerably. Upto a change in Y of about 

in 0 

2.0 Km, the change/is very small. However beyond this point 
any further increase in Y% increases Q in an asymtotic 
manner. At distance of about 20 Km, the 0 values for 
the three different layer heights chosen for simulation 
tend to converge to the same point. The change in Q ■ 
for different heights is 

Tor layer height 2 Km 2,5° to 14.4° 

Tor layer height 5 Km 4.6° to 14*6° 

Tor layer height 8 Km 6.6° to 14.7° 

The increase in 6 reduces the power received by the 
factor ( sin ^ Thus the expected drop in signal 

level is very small foz* Y == 2 Km and pronounced beyond 


this value 
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S. hanges in common vol-ume length ; 

The length of common volnme is maximum at 
the grazing incidence of the antenna heams in the 
mid path plane , Any swing of the heams along and/or 
across the path reduces it. However for a given Y 
position and given layer height, this change is very 
small for along path variations than that for cross 
path variations (Pig. 4 . 5 ), The common volume decreases 
with height. But the change in the length of common 
volume with reference to its spatial position is 
more for lower layer heights than that for higher 
heights. It is seen that for the three layer heights 
chosen for the simulation model the change in length 
of common volume is very small upto about Y=2 jim. 
Purther "beyond about Y=18 Kms the length of common 
volume (and hence common volume) is approximately 
same f on all the height levels considered. Thus 
"beyond this point the reduction in the received signal 
level if any, would not be due to the common volume 
reduction but entirely due to changes in scattering 
angle and such allied factors of the scatter 
integral. 
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The length of the common volume is derived 
^ (5.6) 

from eq_n.^ which is not dependent on scatter plane 
inclination * It can -he said, therefore, that the 
common volume centre cross-section is approximately 
independent of the scatter position within the region 
considered for simulation purposes* • 

Received newer versus the snatial -position ; 

The signal level drops with height and position 
across the signal path (Pig. 4.4 and 4,7 ). The 
apparent width of the layer increases along path and 
decreases across path. 

Por a representative layer of width 400 meters 
and located at a height of about 5.0 Km for exponent 
p = 2, ^ , 5 the variations in received power with 
Y is shown in Pig. 4. 7. The plot reveals that the received 
signal has approximately the same statistical properties 
for p = 2 and In case of p=5, the plot shows appro- 

ximately same statistical properties upto 8 Ems beam 
swinging. However beyond this point, the fall in the 
received power level is slightly more • 
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The value of for p = is more close 
to the observed values from the experimental Kanpur - 
Nainital link. 

The system response for different layers 

thickness is shown in Pig. 4»4. Three layers of 100, 

200 and 400 meters each are considered and signal 

received computed. It is seen that for the layer widths 

100 meter and 200 meters, the signal level is approxi- 

respective 

mately same for ^ heights considered. However 

when the layer thickness is increased to 400 meters, 
at lower height (2 Km), the signal level is more by 
about 1,5 db in the mid path plane and drops down to 
the 100 meter level beyond Y = 7 Kms. This indicates 
that the layer width of 100 meters or so is in a way 
representative of the actual tropospheric common 
volume • 

The simulated signal level in the mid path 
plane for the grazing incidence is -66.9 dEb . fflae 
observed maximum signal level for our link is approxi- 
mately -65 dB'n. Thus there is disparity of about 4 dBn. 
This appears to be due to the fact that the contribution 
due to specular reflection is not included in the 
simulated model . Pur ther in the actual experimental 
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set up large extent of the signal path lies along one 
of the major river. In such a case, the lapse rate 
of moisture with height is expected to he very high. 
This will create favourable conditions for ducting. 
Thus the observed signal is likely to have a component 
due to specular reflection and a component due to 
ducting. This appears to explain the fact that the 
signal strength obtained from simulation is less than 
the ©bseinred value, 

6,2 Conclusion ; 

In this thesis attempts have been made to 
determine the effects of the system geometry and the 
layered atmospheric sturcture of the troposphere on 
the received signal. The simxiLated results 
indicate lower signal strength than the experimental 

results. One of the troposphere parameters function 

2 -12 
0 has been, used with value equal to 5* x 10 
n 

throTjghout the simulation. But this enters linearly 
into the received signal. Therefore any change in the 
signal level by about few d Bo will yield any value 

of signal distribution. 
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For tlie simulation purposes, the medium is 
assumed to he isotropic at the frequency and scattering 
angle used. The shadow-zone of the earth increases 
monotomically in height at greater distance across 
path. Hence extended coverage of upto 20 Kms beam 
swinging used for simulation may not give accurate results 
on the actual link. 

^•5 Recommendations ; 

The simulation model developed in this thesis 
does not indude speculah reflection contrihutions. The 
mathematical treatment given in chapter V gives analytical 
expression for this contribution. 

Also the scattering layers are assumed to be 
horizontal. The effect of the presence of the tilted 
layers within the common volume has not been considered.. 

The atmosphere has been assumed to be isotropic for the 
wavelength and scattering angle used, which may not be true 
in actual case. 

Hence the simulated model could bo further extended 

to take care of 

a) Anisotropy of the medium 

b) Contribution due to specular reflection and 

c) Effect of^presence of one or more tilted layers 
within common volume. 
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/ JOB 
/ FOR HAI 

iOCSiCARD^lYpFWRITER) 

ONE KORU iNTEGERS 

NON FR.iCEJS PROGRAM • 


simulation program for Calculation of the power recieveo from the 

TROPOSPHERE CONSiDERING TROPOSPHERETO BE LAYERED STRUCTURE. 


THE SIMULATICN PROGRAM HAS BEEN RUN ON IBM 1800 COMPUTER. 

THE BASIC ASSUMPTICHS MADE FUR THE PURPCSES OF ANALYSIS ARE 

DTHE EARTH IS ASSUMEO TO BE A SPHERE OF RADIUS 4/3 TIMES 
THE ACTUAL RADIUS TO TAKE CARE OF THE DECREASE IN THE 
VALUE OF REFRACTIVITY WITH INCREASE IN HEIGHT. 

2)THE RECItVER AND TRANSMITTER ANTENNAS ARE ASSUMED TO BE 
LOCATED ON THE SURFACE OF THE EARTH. 




FOLLOWING ARE THE INPUT PARAMETERS 

2D=THE SEPARATION BETWEEN THE TERMINAL STATIONS* ' / ■ r " 
R=EFFECTIVE RADIUS OF THE EARTH WHICH IS TAKEN AS 8500 KMS. , ; 

PT = TRANSMITTED POWER IN WATTS, '' — ■'■—L;'*; 

ALAMO = THE SYSTEM WAVELENGTH BEING USED ■■ '/ 

BeETA=HALF BEAM-WIDTH OF THE ANTENNA BEAM WHICH IS' . 950 DiGRf ES.'L*:,' • 
HT=THE HEIGHT OF ThE SCATTERING LAYER. THREE DIEEERENT HEIGHTS ' 
NAMELY 2,5, AND 8 KM ARE CONSIDERED FOR THE PURPOSES, OF THE SlMULATI' 

ie« **!{<){{:}:*%:{!*!{£*:<! sji# ^ « 3?=^ :{«:?! :j£i Jit s)t5}!*% ##14! :«:3jS*S{« *{«!»« 


PROGRAM STARTS nU * ■ ' ? 

REAL NT,NR, LT, LR,NTMAX,NRMAX,NETMX,NETMNi,NEETA „ 
INTEGER TN ■ . 

DIMENSION PRWATISO), PRDBStSO) ' 

DIMENSION H 1 ( 1 0 3 , H2 ( 1 0 J ■' f 

C OM MON/ / X P , Z ET AP , NE E T A, ZTL0W\ » THE T A , ■ ‘#;0; : 

TN-1 ' ■ " ' ' 

ALAMD=. 14826 
PT=1000. 



R=8.'>t 0 6 
0=16i>.S 03 


PAI IS THE C0MVEF;SI3N FACTOR FOR COf'iVERTiNG RADIANS TO DEGREES 
Pl=4.«ATAri{ 1. ) 

PAl = Pi / Irt'J. 

D^LX-4.t: 03 
' DtLY = ‘300. 

aCETA IS THE HALF BE AN-WIOTH OF THE ANTENNA BEAM. 

6£ETA=PAI*.95/2. 

CM IS THE STRUCTURE FUNCTION COEFFICIENT. 

Crj=S.E-12 
DELZ =100.0 
HT=2.c 03 
H0=0 .S«D*i'D/R 
READ! S, 127) FO 
27 FORMAT (F6. 5) 

NOL DEFINES THE NUMBER OF THE LAYERS BEING CONSIDERED 

■ f ' 

N0L=1 

DO 31 NN=ltNOL 

ATKNS IS THE THICKNESS OF THE LAYER BEING CONSIDERED 
ALOBD IS THE LOWER BOUNDARY OF THE LAYER “ ‘ 

REAO( 5,32) ATKNS 
ALOBC=HT+HO 
Hl( NN) =ALOBD 
H2(NN)=H1(NN)+ATKNS 
31 CONTINUE 

32 FORMAT! 2X, E15. 8) 

XP,YPAND ZP DEFINE TEE ARBITRARY CENTER OF THE COMMON VOLUNE^;' 

T'-''r 

8 XP=0.0 

9 YP = 0.0 
J = 1 

ZP=HPRIM+HT 

C.' ■ : . 

HPRiM IS THE OISTANCF, OF THE EARTH *S ' SURFACE .FROM THE X-Y .PLANE FOR. 

", 'M' F.' -'I 

5 HPRIM=HO-0. 54{ XP^XP+YP*YP)/R ^ ^ ' 

X,Y AND 2 DEFINE ANY ARBITRARY ROIN-I* IS'ASSUMfcO,- ;T0 

COINCIDE WITH THE CENTER 0P’:'’fWE; ' v- 

X=XP 

Y=YP 



no o o o o n o o o n o o o o o o 


1 = IP 

DFN=SQRT{YP*YP+Zf'*ZP ) 

PRWAT(J) IS THE KtClEVED PUWFR IM MILLIWATS 

PUNC IS THE VALUi OF THE INTEGRAND IN THE SCATTER INTEGRAL 

PRWAT( J }=0.0 
FUNC=0.(< 

WRrrE(TN,112) J,X,Y,HPRIM 

112 FnRMAT(2X,3HJ= , 13 ,2X,6HX = , FIS .8 , 2 X , 6HY =,E15.8, 

1 2 HhHPRi M - f E 1 5. 8/ 2X j 3 , ‘>X » , 1 7X » , 1 7X » 6H>>=***^* ) 

WRITE (T.N, 113) 

13 FORMAT (BX, ‘THETA IN' , 1 IX , • G ANA I M • , 1 2 X , ‘PR I N • 1 1 1 X , * PR I N‘ /BX , 

I ‘DEGREES', 12X, 'DEGREf;SM2X, ‘WATTS* ,12X, ‘DBM'/J 
DO 2U0 K=l,40 
IF(YP) 22,22,23 

22 GAMA=89.475*PAI 
GO Tu 24 

23 GAMA=ATAN( ZP/YP}-BEETA 

24 CONTINUE 
GAMA=GAMA+PAl*FLuAT( K) /40. 

GMADG=GAMA/PAI 

ZETAP=YP«C0S(GAMA)+ZP*SIN{GAMA) 

THETR AfWD THETT ARE TFE AGLES MADE BY THE BORES IGHTS OF THE 
RECIEVER AND TRAtlS YI TIER ANTENNA RESPECTIVELY. 

THETT=ATAN{ ZETAP/( C+XP )) 

THETR=ATANUETAP/{ C-XP) ) 

THETA IS THE SCATTERING ANGLE. 

THeTA=(THETT+THETR )/PAI 

THIS PROGRAH CALCULATES THE SLOPES OF THE ANTENNA CONE TANGANT 
PLANES AND HENCE THE X AND ZEETA LIMITS OF THE ONTEGRATION. 

-•T ' 'i > t. ■ ^ ■ 

TRl-SINITHETR+BEETAI /COSI THETR+BEETA) 
TR2=SIN{THETR-8EETA)/COS{THETR-8EETA) 

TTl = SIN(THETT+0eETA) /C0S(THETT+8eF.TA) 

TT2=SIN(THETT-8eeTA) /C0S(TH£TT-36ETA) 


XMIN= D*{TR2-TT1)/ (TR2+TTL) 

;2 XMAX= D«‘(TRi-TT2)/ {Tra+TT2) 
ZTMIN= 2.*D*TR2*1T2/ {TR2+''T2) 
ZTMAX= 2.*D*TR1>!'TT 1/ ( TRl+TTl ) 



oooo non ooooooooo ono 


THIS program ENSURcS THAT THE INTEGRATION STARTS ONLY WHEN jHt 
SCATTER PLANE COHES WITHIN THE ANTENNA CONES COMMON ZONE. 

IF( YP) 26,26,27 

26 EXPR==ZTMAX^CaS {GAMA)-ZP*SIN( aeETA) /CQSI BEETAI 
GO TO 28 

27 EXPR = ZTMAX«C0S {Gam A) -(YP + ( ZTMAX-ZTMIN )/2. ) 

28 ir-(-XPR) 7,2C0,2i:Q 
/ CriNTl.NU.- 

ALNTH IS THE LENGTH DF THE COMMON VOLLNE 

ALNTh=XMAX-XMI N ' 

ZEETA=Y*COS( GAMA )+Z* SI N( GAMA) 

THIS program CALCULATES THE VALUE OF THE POLARIZATION FACTOR 
FUR KANPUR NAINIiAL LINK THE A-)TENNA HORN FEEDS ARE VERTICALLY 
POLARIZED. HENCE THE EXPRESS ION CALCULA TED FOR ACCOUNTING OF 
POLARIZATION EFFECTS IS FOR VERTICAL POLARIZATION . 

AK1 = Z=S‘Z*{ ( D + X) =!>«2+Y«Y3 
AK2=( {D+X)*=>=2 + Y=i=Y) «*2 
AK3= { D-X) «X=2 +Y *Y -kZ ^Z 
AK=4.’S=Z=«'Z*D=«=D*{D+X )**2 
SKSa=l.-AK/{ {AK1+AK2 )*AK3) 

LT' AND LR ARE THE SQUARED DISTANCES OF POINT P FROM TX AND RX. 

■ ■„ . . ;; » . ■ i 

LT=(D-»-X) **2+ZEETA*ZEETA 
LR*(D-X)*( D-X) +ZeETA>!‘ZeETA 

THIS PROGRAM DETERMINES THE NEETA LIMITS OF THE INTEGRATION 
THE NEETA LIMITS ARE CALCULATED FOR THE VALUES OF NEETA GREATER 

NT=8EETA*0EETA"I ATANCZEETAZID+X) }-THETT)**2 
IF(NT“0.03 200,200,1 

1 NTMAX=SQRTCNT*LT3 

NR=8EETAX‘BEETA~I ATAN C 2E6TA/{D~X3 }-THE IR )X=#2 
IF{NR-OQO) 2OO»20O,2 

2 NRMAX=SQRTINR*LR) 

IF! NTMAX-NRMAX ) 15,26,25 
15 NETMX^NTMAX 
N£TMN=-N£TMX 
:G0 -TO 35 
'25' 'NETMX=NRHAX 
' : N6THN=-NcTMX 


C THIS PROGRAM SOLVES TPE SCATTER INTEGRAL 

c F!)r scatter integral solving, SIMPSON *n method cf integration is 

c 

35 JJ= 20 ■ 

RJ=FlOAT{ JJ) 

AP=( 2TMAX-ZTMIN) /R J 

ztlow=ztmin 

D>J 5u I =1, JJ 
Z1^ZP-NETMX«C0S( GAMA ) 

Z2 = /-P+NETMX*C3S( GAf^A ) 

DO 45 LL=1,N0L 
IF(HI(LL)-Z2) 61,45,45 

61 IF(H2{LL)-Z1) 45,45,62 

62 IF(H1(LL}-Z1I 64,67, 67 
64 ALLT=Z1 

GO TO 6S 

67 ALLT=H1{ LL) 

68 IF{HZ( LL)-Z2) 36,36,37 

36 UPt.T=H2{LL) 

GO TO 29 

37 UPLT=Z2 

29 W=UPLT-ALLT 

AR=(NETMX-NETMN) /RJ 
NEETA=NETMN 
DO 45 KK=1, JJ 

IF (K-I) 16,16,17 

17 YNU=ZETAP*COS( GAMA ) + NeTMN*SIN{ GAMA } 

IF(YLMT-YNU) 45,16,16 
16 CONTINUE 

Z = ( ZTLOW’S'ZP-NEETA^YP ) /DEN 

C HS IS THe;-,HEIGHT OF TH,*e SCATTERING LAYER 

c ■ .V - 

HS=i-HPRiw 
delh=z~allt ' 

I F C Q§tH^0*:0 1 ' 

6 IFIOiLH~WJ 
3 CONTINUE -t ^ ^ 

h^^ixmax-xhinI/rj 

XI^XHIN ’ 

NUM*4J/2T.-. .. 

=M=1#NUH 

#ilNC«FUNC+H*( F I X n +4 .X'FI Xl+H) + F{ XI +2. H /3. 

PftWATI J) = PRWAT( J)+l . 91^.000 L’!'CH*PT--!'SKSU*FUNC*ALAMD«*{5,/3. )^IQ)Z* 
|F(PRWAT(J)) 200,200,21 
21 PRDBSI J ) = 10,’i=AL0G( PRwAT( jn / ALOGI 10. ! 

40 XI=XI+2.*H 


45 HCET A=NEETA+AR 
50 ZTLGW=ZTLOW+AP 

YLf^T-ZETAP^COS (GAhA) +NETMNi=!'SIN( GAMA} 

200 CONTI'AJE 

C , • 

c 

c 

WRIT>(T'^,1L4) THcTA,GFADG,PRWAT(J) ,PR0BS{JI 
114 Fi.lRMAK 5X,4(E1 5. 3X ) ) 

WRIT!: (TN, 129 ) ALf Th , THETA 

129 F0RMAT{ 2X, 'LEMGTH OF COMMON VOLUME = • , E 1 5. 8 , 2X , • THE TA = * , E 1 5. 8/ ) 
IF(YP-20.E 03) 55, 56 ,56 
55 Y=y+0£LY 
GO TO 5 

56 IF{XP-20.E 03) 57, 58,58 

57 XP=XP+DELX 
GO TO 9 

58 IF{HT-8.E 03) 59,60, 60 

59 HT=HT+3.t 03 
GO TO 8 

60 STOP 
END 

// FOR 

FUNCTION F{ X ) 

REAL NEETA,LTtLR 
INTEGER TN 

C0MM0N//XP,2ETAP ,NEETA,ZTLOW , THETA FC 
C 

c 

TN=1 

PIc 4.:«(ATAN( 1. ) 

0=165. E 03 

THETT=ATAN( IETAP/I D«-XP| ) 

THETR=ATANI ZET W/.fD-XPl ) 

T8=NeETA*N6ETA/C 1 **2^271 OW’i'ZTLOW ) 

TB8={ ATANf ZTti3W/ID#X ) J-TrtETT )**2 
BetAT=sS.QRT«T8+T86) 

C vf. , .. 

C ' 

TR^||IE8t'A*NEETA/{ JD-X )**2+ZTLOW=«'ZTLOW ) 

I AT ANT Z T LO W/ ( 0 -X ) ) -T HE TR ) * « 2 
;#'fAR=SQRT{TR+TRR) 

C.’ 

C', 


PA1=P1/180 



B£tOG=B£TAT/PAI' 

'88 B£R0G=8ETAR/PAI 
c 

C THEMAXIMUM value of the ANTEflNA GAIN HAS SEEN USED AS 42 DBS. 

GTGR= LXP(-F 0 *« BeTDG=<'BETDG+BeRDG«BERDG n (8.4 ) 

LT=( D+X) >!'*2 + ZTLOW>i'*2 
LR=(iJ-X)=i=*2+ZTLOV.=i=^2 
C 

ThTRD=ATAN( ZTLOW /( !>X ) ) + ATAN( ZTL3W/( D+X H 

C 

A = SI (THTRO/2. )**( 11 , /B. ) 

C 

c 

F=GTGR/(LT*LR*A) 

RtTU ['■!'( 

END 

// XEQ MAIN 
*CCEND 
$ENTRY 
C 

G THIt ABUV6 PROGRAM (JHICH HAS BEEN DEVELOPED FOR IBM 1800 MACHINE 

e CAN VERY CONVENIFNTLY CONVERTED TO BE USED ON IBM 7044 BY THE 

C CHANGE OF THE COl.TROL CARDS AND THE FCRMAT STATEMENTS SUITABLY. 

C. ^ ■ ■ • 

C use SUITABLE CONTROL CAROS. 

C USE SUITABLE DATA 

C 
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